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Inelastic neutron scattering is employed to study the reciprocal-space structure and dispersion of magnetic 
excitations in the normal and superconducting states of single-crystalline RbogFe! 6 Se 2 . We show that the 
1 _H recently discovered magnetic resonant mode in this compound has a quasi-two-dimensional character, similar 

i— H to overdoped iron-pnictide superconductors. Moreover, it has a rich in-plane structure that is dominated by 

four elliptical peaks, symmetrically surrounding the Brillouin zone corner, without a/5 x a/5 reconstruction. 
£N| We also present evidence for the dispersion of the resonance peak, as its position in momentum space 

depends on energy. Comparison of our findings with the results of band structure calculations provides 
strong support for the itinerant origin of the observed signal. It can be traced back to the nesting of 
electronlike Fermi pockets in the doped metallic phase of the sample in the absence of iron-vacancy ordering. 

PACS numbers: 74.70.Xa, 74.25.Ha, 78.70.Nx, 74.20.Rp 
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The newly discovered iron selenide superconductors 
A x ¥e 2 - y Se 2 (A = K, Rb, Cs) [1] became famous for their rel- 
atively high critical temperature, T c = 32 K, observed con- 
currently with a strong antiferromagnetic (AFM) order that 
persists far above room temperature [2, 3]. However, an or- 
dered moment as large as 3.3 jU B /Fe [4] renders microscopic 
coexistence [5] doubtful. The superconducting (SC) phase 
usually appears in samples close to the 2:4:5 stoichiometry 
[3, 6], which is at the same time the optimal composition 
for the ordering of Fe vacancies into a a/5 x a/5 superstruc- 
ture, grouping the occupied iron sites in plaquettes of four 
ferromagnetically aligned moments. On the one hand, ex- 
periments [3, 7, 8] and band structure calculations [9, 10] 
suggest this superstructure to be insulating. On the other 
hand, angle-resolved photoelectron spectroscopy (ARPES) 
revealed a Fermi surface (FS) dominated by a large electron 
pocket at the M point [11, 12]. Recent reports reconcile 
these seemingly contradictory findings by the observation 
of several coexisting phases, seen in transmission electron 
microscopy (TEM) [13-16], scanning tunneling microscopy 
(STM) [7], x-ray diffraction [17], ARPES [18], magnetiza- 
tion measurements [19], muon-spin rotation (/iSR) [20], 
Mossbauer [21] and optical [8] spectroscopies. While STM 
studies observed the SC gap on a vacancy-free surface [7], 
TEM measurements suggested that the second phase is an 
iron-vacancy disordered state [13, 14]. The phase separa- 
tion scenario clearly needs more clarification in terms of the 
structure and stoichiometry of the SC phase for a consistent 
understanding of these observations. 

In iron pnictides, it is established that the SC order pa- 
rameter changes its sign between the hole- and electronlike 
sheets of the FS [22-24]. Despite the absence or strong 
reduction of the hole Fermi pocket in iron selenides [25], 
different kinds of a sign-changing gap have also been sug- 
gested [26-28]. The recent finding of a magnetic reso- 
nant mode in the low-energy spin-excitation spectrum of 
Rb 8 ¥e 1 6 Se 2 below T c [29] supports these unconventional 
pairing scenarios. Its wave vector, Q = (\\\), can be 



reconciled with theoretical calculations performed for the 
electron-doped phase with the A x ¥e 2 Se 2 stoichiometry (not 
matching the average chemical composition of the sample) 
under the assumption of a d-wave symmetry of the SC order 
parameter [27, 28]. Alternatively, the metallic phase could 
possibly be associated with (i) a vacancy-disordered struc- 
ture [13, 30], which would effectively result in a rigid-band 
shift and broadening of the A^Fe 2 Se 2 electronic bands [10], 
(ii) an electron-doped A x Fe 4 Se 5 phase [31] with full or par- 
tial vacancy ordering, or (iii) possess a different a/2 x a/2 
superstructure that corresponds to the A x Fe x 5 Se 2 composi- 
tion [13, 16, 32]. However, to the best of our knowledge, 
first-principles calculations of the spin-excitation spectrum 
are not yet available for any of these alternative scenarios. 

To be able to differentiate between the mentioned possi- 
bilities and thus try to verify the origin of the spin-excitation 
spectrum, we have performed a detailed study of the 
reciprocal-space structure and the dispersion of the pre- 
viously reported resonant mode. We show that the resonant 
magnetic excitations in the SC state of Rb x Fe 2 _ y Se 2 are 
quasi-two-dimensional (2D) and exhibit a complex in-plane 
pattern, dominated by four elliptical peaks that symmetri- 
cally surround the corner of the unfolded Brillouin zone 
(BZ) [33]. This result is consistent with the dynamic spin 
susceptibility of an electron-doped A x Fe 2 Se 2 compound, cal- 
culated in the SC state from a tight-binding model of the 
band structure by means of the random phase approxima- 
tion (RPA) [27]. 

The sample for this study is identical to the one used in 
Ref. 29. It comprises several coaligned single crystals with 
the chemical composition Rbg gFe! 6 Se 2 and a total mass 
of ~ 1 g. These crystals are bulk superconductors with a 
T c of 32 K [29], which have been characterized by trans- 
port, magnetometry and specific heat measurements (batch 
BR16 in Ref. 6). The experiments were conducted at the 
thermal-neutron spectrometer IN8 (ILL, Grenoble), which 
was operated both in the triple-axis-spectrometer (TAS) and 
in the FlatCone multianalyzer configurations. The latter al- 
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Fig. 1 (color online), (a) Longitudinal momentum profiles through 
Q = {\ \ — |)in the normal state at 35 K and in the SC state at 
1.5 K, both at E = 14meV (b) L-dependence of the intensity at 
E = 14meV along Q = (^ \L) at the same temperatures, (c) The 
respective difference of the two signals fitted to the Fe 2+ magnetic 
form factor (solid line), (d) Energy scans at the resonance position, 
measured at L = —0.5 with k { = 4.1 A -1 and at L = 0.5 with 
k f = 2.662 A -1 in the normal and in the SC states, (e) Difference 
of the SC- and normal-state intensities at integer and half-integer 
L. The solid line through the data points for k { = 4.1 A -1 is a fit 
with a Gaussian function superposed on the difference of the Bose 
factors for both temperatures (dashed line). The data points taken 
with k f = 2.662 A -1 are shifted down by 150 counts for clarity. 



lowed us convenient mapping of the reciprocal space at a 
constant energy. In the FlatCone configuration we utilized 
Si(lll) monochromator and analyzer with the fixed final 
wave vector k f = 3 A -1 . For the TAS measurements, we 
used pyrolytic graphite (002) monochromator and analyzer 
with double focusing. The TAS measurements were done 
with constant fc f = 2.662 A -1 and 4.1 A -1 , and a pyrolytic 
graphite filter was installed between the sample and the an- 
alyzer to suppress higher-harmonic contamination. In order 
to measure the dispersion and the intensity distribution of 
the resonant mode along the c axis, we mounted the sample 
in the {2HHL) scattering plane. Subsequent investigations 
of in-plane excitations were done in the (HKO) plane. Here 
and throughout the paper our notation is given in unfolded 
reciprocal lattice units (r.l.u.), which refer to the iron sub- 
lattice with the lattice parameters a = b = 2.76 A and 
c = 7.25 A [29]. The existence of an Fe-vacancy-ordered 
AFM phase in our sample has been verified by measuring 
the magnetic superstructure reflections [29] . 

In Fig. 1 (a) we present longitudinal [as seen from (| | 0)] 
momentum scans at the resonance energy, Hco res = 14meV 
along the (| K — ~) direction in the normal state at 35 K 
and in the SC state at 1.5 K. Already in the normal state, 



we observe a substantial magnetic response, which becomes 
considerably enhanced below T c . The center of both peaks 
lies at K = (0.244 ± 0.002), close to the commensurate po- 
sition at K = \. In Fig. 1 (b), the INS intensity atQ = (| \) 
and E = 14meV is plotted as a function of the out-of-plane 
momentum component along (| ^ L) for both temperatures. 
The intensity difference between 1.5 and 35 K [Fig. 1(c)], 
representing the resonant enhancement, is maximized at 
L = and then monotonically decreases for larger L follow- 
ing the Fe 2+ magnetic form factor. This closely resembles 
the 2D nature of the signal in overdoped BaFe 2 _ x Co x As 2 
[23], but is in contrast to the behavior of underdoped 
BaFe 2 _ x Ni x As 2 , where it is modulated as a function of L 
and exhibits a maximum at L = | [33]. 

The weak L-dependence of the resonant signal can also 
be inferred from the energy scans measured above and 
below T c , such as those presented in Fig. 1 (d) for fc f = 
2.662 A -1 and 4.1 A -1 . The intensity difference between the 
two temperatures is shown in panel (e) for both integer and 
half-integer L . Neither the energy nor the amplitude of the 
resonance peak shows any notable L-dependence beyond 
the uncertainty of the fit. This fact is consistent with ARPES 
measurements of the weak k z -dispersion of the electron 
band at the M point and of the SC gap [12]. 

Based on the quasi-2D character of the magnetic intensity, 
we have mapped out the resonant enhancement of spin 
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Fig. 2 (color online), (a) Color map of the reciprocal space, show- 
ing intensity difference between the SC and normal states at 
E = 15meY measured by the FlatCone detector, (b) The same 
map as in (a), rebinned on a 81 x 81 grid, symmetrized with re- 
spect to the mirror planes and smoothed using a Gaussian filter 
with 1 pixel standard deviation, (c) Longitudinal cuts (along the 
short axis of the ellipse) through the data in (a) at all four reso- 
nance positions. The intensity was integrated within a window of 
0.28 r.l.u. in the direction perpendicular to the cut. (d) The same 
in the transverse direction (long axis of the ellipse) . The intensity 
was integrated within a window of 0.12 r.l.u. in the direction per- 
pendicular to the cut. The widths of the integration windows are 
given by the horizontal and vertical bars in panel (a), respectively. 
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Fig. 3 (color online). TAS-mode data, (a) Longitudinal momentum 
scans through the center of the ellipse at Q x = ( | \ 0) (triangles) 

and Q 2 = ( \ | 0) (squares and diamonds) as indicated in sketch 
(c) at E = 15 meV The intensity in the SC and in the normal states 
(top) is shown together with their difference (bottom) . (b) The 
same for transverse momentum scans at Q 2 = (| \ 0). (c) A frag- 
ment of the FlatCone map from Fig. 2 that illustrates the directions 
of the scans shown in this figure. Panels (d) and (e) show only 
the difference in intensity between SC and normal states, (d) Mo- 
mentum scan at E = 15 meV parallel to the longitudinal direction 
at K = 0.42, offset from the center of the ellipse, (e) Momentum 
scans at different energies along the short axis of the ellipse. The 
plot at E = 15 meV is an average of the two profiles in panel (a) at 
both resonance positions. 



excitations at E = 15 meV in the (HKO) scattering plane by 
means of the FlatCone multianalyzer. Figure 2 (a) shows the 
difference of intensity maps measured around the BZ corner 
in the SC and normal states. We observe resonant intensity 
at all four symmetric positions equivalent to (~ - 0). In order 
to reduce the statistical noise in the data, we have rebinned 
this data set on an 81x81 grid and symmetrized it with 
respect to four mirror planes of the reciprocal space, with 
subsequent Gaussian smoothing. The resulting intensity 
map is shown in Fig. 2 (b) as a contour plot. 

One sees that the in-plane shape of the resonant intensity 
takes an elliptical form, elongated transversely with respect 
to the vector connecting it to (| 1 0). We emphasize this by 
presenting cuts through all ellipses in the map of Fig. 2 (a) 
in the longitudinal [Fig. 2 (c)] and transverse [Fig. 2 (d)] 
directions. The intensity is integrated over the whole exten- 
sion of the ellipse perpendicular to the cut, as indicated by 
the black bars in (a), in order to capture the whole resonant 
intensity. We observe an agreement between equivalent cuts 
with the same orientation. The ratio of the peak widths 
in the transverse and longitudinal directions results in an 
aspect ratio of 2.1 for the resonance feature. 

Next, we present momentum scans measured by TAS in 
longitudinal [Fig. 3(a)] and transverse [Fig. 3(b)] direc- 



tions through the ellipse in the SC and normal states. Again, 
strongly anisotropic widths of the transverse and longitu- 
dinal profiles are observed in the intensity difference. The 
peak in the longitudinal direction for both resonances near 
(| \ 0) and (| | 0) in Fig. 3 (a) is found at an incommen- 
surate position of H = 0.78 or K = 0.78, respectively, as 
marked by the arrow. This is also consistent with the Flat- 
Cone data in Fig. 3 (c), where the peak intensity is offset to 
the right from H = |. 

An elliptical in-plane shape of the resonance has also been 
observed in BaFe 2 _ x Co x As 2 [24, 33] and in Ba 1 _ x K x Fe 2 As 2 
[34] at the BZ boundary, so that both axes of the ellipse are 
aligned along the natural mirror planes of the reciprocal 
space. In Rb x Fe 2 _ y Se 2 , however, the ellipse could be asym- 
metric, because H = | is not a natural high-symmetry plane. 
Indeed, the shape in Fig. 3 (c) suggests a slight bending of 
the ellipse towards (| \ 0). In the colormap in Fig. 3 (c), 
we also observe weak streaks of intensity reaching towards 
(| 1 0) and (| | 0), barely above the statistical noise level, 
which could form parts of a ring connecting all four reso- 
nance positions. Nevertheless, the peak profile measured 
parallel to the longitudinal direction and offset by 0.08 r.l.u. 
from the center of the ellipse [Fig. 3 (d)] does not show 
any notable shift of the peak center beyond statistical un- 
certainty. This indicates a nearly symmetric (non-curved) 
shape of the resonance peak in the vicinity of its maximum. 

Finally, we turn to the in-plane dispersion of the reso- 
nance, which could be studied due to the broad distribu- 
tion of the resonant intensity in energy, as can be seen in 
Fig. 1 (e) . Figure 3 (e) presents longitudinal momentum 
scans of the resonant intensity at 12, 15, and 18meV Here, 
the peak center shifts from H = (0.764 ± 0.002) r.l.u. at 
12meV to H = (0.782 ± 0.003) r.l.u. at 15meY although 
we do not resolve a further shift upon changing the energy 
to 18meV Moreover, comparison of the peak position at 
L = -0.5 [Fig. 1 (a)], centered at H = 0.244 ± 0.002, and 
at L = [Fig. 3 (a)], where it is shifted to a position equiv- 
alent to H = 0.218 ± 0.003, also suggests a small (~ 10%) 
variation in the peak position along the c axis. 

To verify the origin of the observed spectrum of spin 
excitations in Rb x Fe 2 _ y Se 2 , we will now compare our ex- 
perimental observations with the results of band structure 
calculations. For this purpose, we employ the tight-binding 
model that was introduced in Ref. 27 to describe the elec- 
tronic structure of an electron-doped A x Fe 2 Se 2 . The chemi- 
cal potential has been adjusted by a rigid-band shift of the 
bands to match the positions of the magnetic resonant peaks 
in the calculated susceptibility with the experimental data. 
This resulted in a doping level of ~ 0.18 electrons/Fe. To en- 
able direct comparison between the theory and experiment, 
we have calculated the imaginary part of the dynamical spin 
susceptibility at the resonance energy, #"(Q, co res ), both for 
the SC and the normal states, as described in Ref. 27. For 
the calculation in the SC state, we have assumed a d x 2_ y 2 
gap A(k) = A (cosk x - cosk y ) [35]. The color map in 
Fig. 4 (a) shows the respective difference of the two quan- 
tities, ^s C (Q, co res ) - ^(Q, co res ), within the (HKO) plane, 
isotropically broadened by a Gaussian resolution function 
with a standard deviation of 0.02 r.l.u. Comparison with the 
experimentally measured resonant intensity map in Fig. 2 (b) 
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Fig. 4 (color online), (a) The difference of the calculated imaginary 
parts of the dynamic spin susceptibility for the SC and normal 
states, taken at the resonance energy #g C (Q, ^> res ) — Xn(Q» ^res)- 
The calculation was done within RPA from the tight-binding band 
model of A x Fe 2 Se 2 [27], which was rigidly shifted to match the 
experimental peak positions. An isotropic Gaussian broadening 
with a standard deviation of 0.02 r.l.u. has been applied to mimic 
the experimental resolution, (b) The resulting FS in the (H K 0) 
plane corresponds to the doping level of 0.18 electrons/Fe. The 
black arrows are the in-plane nesting vectors responsible for the 
resonance peaks observed in our present study. 

reveals good agreement between the two Q-space patterns, 
as both the orientation and the aspect ratio of the elliptical 
peaks is well captured by the calculation. The origin of these 
peaks can be traced back to the nesting of electronlike Fermi 
pockets, as indicated in Fig. 4 (b) by black arrows. 

To conclude, the fact that the complicated pattern of reso- 
nant intensity in Q-space can be successfully reproduced by 
our calculation strongly supports the itinerant origin of the 
observed magnetic response. The signal shows no signatures 
of the a/5 x a/5 reconstruction, indicating that it originates in 
the metallic phase of the sample without iron-vacancy order- 
ing, as suggested recently [13, 14, 36]. This distinguishes 
the newly observed signal from the previously reported spin- 
wave excitations in this class of compounds [37] that stem 
from the magnetic superstructure Bragg positions in the in- 
sulating vacancy-ordered phase and are insensitive to the SC 
transition. The incommensurability of the resonance peak, 
as well as its variation with the out-of-plane momentum 
component and with energy, further indicates that it is not 
pinned to a particular position in Q-space, but is arbitrarily 
determined by the level of electron doping, in line with the 
assumptions of Ref. 27. 

Furthermore, we note that iron-pnictide compounds gen- 
erally exhibit a tendency towards a 2D behavior of spin 
fluctuations with an increase of the doping level or T c 
[33, 38]. For example, optimally doped Ba 1 _ x K x Fe 2 As 2 , 
which has the highest known critical temperature among 
all 122-compounds, shows almost no dispersion of the reso- 
nant energy, Hco res , along the c direction [34]. Our data on 
Rb x Fe 2 _ y Se 2 with a comparable transition temperature are 
fully consistent with this trend. 
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